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We investigate controlled phase separation of a binary Bose-Einstein condensate (BEC) in the 
proximity of mixed-spin-channel Feshbach resonance in the \F = 1, ttif = +1) and \F = 2, ttif = — 1) 
states of Rb at a magnetic field of 9.10 G. Phase separation occurs on the lower magnetic-field 
side of the Feshbach resonance while the two components overlap on the higher magnetic-field side. 
The Feshbach resonance curve of the scattering length is obtained from the shape of the atomic 
cloud by comparison with the numerical analysis of coupled Gross-Pitaevskii equations. 
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I. INTRODUCTION 

Ultracold atomic gases provide an attractive testing 
ground for studying dynamics of multicomponent quan- 
tum fluids. It has been shown that dual-species quantum 
gases PHI, two-component Bose-Einstein condensates 
(BECs) comprised of two different hyperfine states [4 y , 
and spinor BECs with different Zeeman sublevels [6|-l8| 
exhibit a rich variety of dynamics. The controllability of 
the intra and intercomponent interactions via Feshbach 
resonance [9[ creates numerous possibilities and enriches 
the physics of multicomponent quantum fluids. 

Miscibility between different components is crucially 
important to the dynamics of multicomponent systems. 
Phase separation in immiscible two-component BECs has 
been studied in Ref. [10]. Immiscible two-component 
BECs have been predicted to have interface instabili- 
ties: the Kelvin-Helmholtz instability in the presence of 
shear flow [11[ and the Rayleigh-Taylor instability jl2| . 
On the other hand, the miscibility between different spin 
components plays a key role in coherent spin dynamics, 
such as the Josephson oscillation [13|, LL4[, spin echoes 
J15L the Ramsey interferometer [5|, spin entanglement 
[161 LL7j , and determination of the magnetic ground state 
of a spinor BEC . 

Papp et at. [2j recently tuned the miscibility in binary 
BECs of 85 Rb- 87 Rb controlling the intracomponent in- 
teraction of 85 Rb. The present system differs from that 
in Ref. [2j in that the miscibility of different spin states 
of an identical species is controlled by intercomponent 
Feshbach resonance. For a binary mixture of two differ- 
ent internal states, the populations of the components 
can be altered at any stage in an experiment via spin 
manipulations. Although the Feshbach resonance has 
been reported on different internal states of an identi- 
cal species jl9J-[21j, control of their miscibility has not 
been discussed. 

The scattering leng th determines the properties of ul- 
tracold collisions [22]. Spectroscopic methods for deter- 
mination of the scattering length by observation of en- 



ergy shift have been demonstrated [23-26]. While these 
methods have advantages in accuracy, they are applicable 
only to the states between which spectroscopic transition 
is available. In addition, the spectroscopic methods have 
a disadvantage for high density system, since an energy 
shift caused by atomic density lowers the precision of the 
estimation. 

In this paper, we control phase separation via Feshbach 
resonance between internal spin states, \F = l^mp = 
1) = |1) and \F = 2,ra F = -1) = |2), of 87 Rb. The mis- 
cibility of these two components is found to depend sensi- 
tively on the strength of the applied magnetic field near 
the Feshbach resonance. The experimental results are 
compared with numerical simulations of coupled Gross- 
Pitaevskii (GP) equations. The excellent agreement be- 
tween the experimental results and the numerical sim- 
ulations allows us to estimate the scattering length be- 
tween the internal states. The experiments and simula- 
tions for various values of magnetic field yield the reso- 
nance curve of the scattering length near the Feshbach 
resonance. Thus observation of phase separation dynam- 
ics can be used as a new method to estimate scattering 
lengths of multicomponent BECs, which does not rely 
on spectroscopic transition and can be used in the high 
density regime. 

This paper is structured as follows. In Sec Jill we in- 
troduce the mean- field formalism for a binary BEC. In 
Sec lIII( our experimental apparatus and conditions are 
described. In SecHV} experimental results are compared 
with numerical simulations. Sec El is devoted to the con- 
clusions. 



II. TWO COMPONENT CONDENSATES 

The dynamics of a binary BEC with inelastic two-body 
losses is described by coupled GP equations [4], [H, [l8j , 
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where ipi is the macroscopic wave function for the \i) state, 
m is the mass of 87 Rb, and V is the trap potential. The 
interaction coefficient is given by g^ = gij—ifiKij/2 with 
gij = 47rh 2 aij/m, where aij is the scattering length and 
Kij is the two-body inelastic collision rate between the 
\i) and \j) states of condensates. The scattering lengths 
of 87 Rb have almost the same values: <222 = 95.00 ae, 
an = 100.4 ae, and ai2 = 97.66 ae = &bg? where <2b is 
the Bohr radius [4| . The scattering lengths determine the 
miscibility of the binary BEC. The phase separation con- 
dition in a uniform system is given by a\ 2 > ttn«22- The 
atom density in two component condensates decreases by 
the two-body inelastic collisions as 
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The two-body inelastic collision rates in the F = 2 
manifold were measured in Ref. [18j, giving K22 = 
1.04 x 10 _13 cm 3 /s. The two-body inelastic loss in the 
F = 1 manifold is negligible and we assume Kn = 0. 

Near the magnetic Feshbach resonance, the inter com- 
ponent interaction gi2 is changed. The interspecies scat- 
tering length in a complex form can be expressed as a 
Lorentzian function [27| . Since the imaginary part of the 
scattering length is incorporated in if 12, we use the ef- 
fective scattering length between the |1) and |2) states 
defined by 
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where the parameters Bo, AB, and 7b are determined 
later. 



III. EXPERIMENT 

Our experimental apparatus and procedure used to 
create 87 Rb condensates are the same as those described 
in Refs. [14j, LL8| except for irradiation of a two-photon 
7r/2 pulse between the hyperfine states. A BEC contain- 
ing 10 6 atoms in the |2, 2) state is created by evaporative 
cooling with frequency sweeping of an rf field for 18 s in a 
magnetic trap. The BEC is loaded into a crossed far-off 
resonant optical trap (FORT) at a wavelength of 850 nm 
and 3 x 10 5 atoms remain in the FORT. The potential 
depth of the crossed FORT is estimated to be about 1 
jiK and the radial (axial) trap frequency is measured to 
be 141 Hz (21 Hz) using the parametric resonance. Af- 
ter sudden inversion of the quantization axis, the |2, — 2) 
state is transferred to the |2, — 1) = |2) state by inducing 
the Landau-Zenner transition using rf irradiation with 
an external magnetic field of 20.5 G. Half the atoms in 
the |2) state is then transferred to the |1,1) = |1) state 
by irradiation of a tt/2 pulse of rf and microwave field 
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FIG. 1: (Color online) (a) Magnetic- field dependence of the 
atom number for evolution times t ev of 25, 50, and 75 ms. 
(b) Inelastic collision rate estimated by Eq.(|2]) at t ev = 25 
ms. Solid lines are Lorentzian functions fit to the data with 
center at 9.10 G. 



for 5 ms. After time evolution for t ev in a precisely con- 
trolled magnetic field B ev , the crossed FORT is turned 
off. The magnetic field is kept on for the first 5 ms of the 
time-of-flight (TOF) to maintain the scattering length 
during expansion of the atomic cloud. After applying a 
Stern-Gerlach pulse, absorption images are obtained. In 
order to take the image of each component, we apply a 
sequence of two short pulses after a TOF times of 15 ms 
for F = 2 and 18 ms for F = 1 with repumping to F — 2. 
The fluctuation in the number of atoms is estimated to 
be 10%. The relative population fluctuates between 0.45 
and 0.55 in each run of the experiment |18l l28j. 

The magnetic-field strength of the Feshbach resonance 
can be determined by measuring the atomic losses. In 
our experimental conditions, atomic losses by one- and 
three-body inelastic collisions are negligible compared 
with that by two-body inelastic collisions [18[ . The total 
number of atoms in the |2) state decreases rapidly with 
increasing trap time t ev because of hyperfine-changing 
inelastic collisions such as |2, — 1), |2, — 1) — >> |l,mi?), |1 
or 2,ra F ) and |2, — 1), |1,— 1) — » \l,m F ), |l,ra F ). For 
a magnetic field far from the Feshbach resonance, the 
two-body inelastic collision rate between the |2) and |1) 
states is estimated to be K12 = 0.5 x 10 _13 cm 3 /s, which 
is obtained by comparing the atom loss in our experi- 
ment and a solution of the GP equation (1). Both the 
elastic scattering length and the inelastic collision loss 
rate are altered in the vicinity of the Feshbach resonance 
[9[ . Figure QJa) shows the total number of atoms after 
the TOF for evolution times t ev of 25, 50, and 75 ms. 
For t ev = 25 ms, density profiles were almost the same at 
both lower and higher magnetic field near the Feshbach 
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FIG. 2: (Color online) Column densities of the |2) and |1) states obtained by the experiment (upper) and the numerical 
simulation (lower) as a function of magnetic field 5 ev . The field of the view is 316 (vertical) x 290 /xm (horizontal). The two- 
body inelastic collision rates if 12 (10 -13 cm 3 /s) and the change in the scattering lengths Aa used in the numerical simulation 
are indicated at the bottom. 



resonance. We estimated Kyi for t ev = 25 ms by solving 
Eq. (J2]) numerically with the single-mode approximation. 
The uncertainty arises from the fluctuation in the initial 
number of atoms. The data are fitted by a Lorentzian 
function as shown in Fig.[TJb), which gives Feshbach res- 
onance field Bo to be 9.100 G and 7b to be 30 mG. The 
fluctuation in the magnetic field B ev is estimated to be 
less than 5 mG by observing the magnetic dipole tran- 
sitions, and the residual gradient magnetic field is esti- 
mated to be 30 mG/cm [14J|. The magnetic-field strength 
is calibrated by Rabi spectra between clock states with 
the microwave irradiation. The uncertainty in this cali- 
bration is 5 mG. The resonant magnetic field obtained in 
our experiment, 9.100(5) G, agrees with the theoretical 
prediction [29j within the experimental uncertainty. 



IV. RESULTS AND DISCUSSION 
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FIG. 3: (Color online) Numerically obtained column density 
distributions for various scattering lengths at an evolution 
time of 75 ms with K 12 = 3 x 10~ 13 cm 3 /s. The relative 
scattering length Aa in unit of as is indicated in the top row. 



The column densities of the binary BEC around the 
Feshbach resonance for t ev = 75 ms are shown in Fig. [2j 
At a magnetic field far from the Feshbach resonance 
(B ev = 8.30 and 10.05 G), the two components exhibit 
phase separation [I, [5[ . The domain structure of phase 
separation depends not only on the scattering lengths 
but also on the number of atoms and the relative pop- 
ulations. The scattering lengths and inelastic collision 
rates for B ev = 8.30 G and those for B ev = 10.05 G are 
almost the same, since their density patterns are similar. 
The behaviors of the binary BEC in the vicinity of the 
Feshbach resonance at 9.10 G change dramatically in Fig. 
[2] due to the change in the scattering length and the in- 
elastic collision rate. The domain structures on the lower 
magnetic-field side near the Feshbach resonance (B ev = 
9.05 and 9.08 G) are quite different from that at B ev = 
8.30 G. On the other hand, the domain structure dis- 



appears at higher magnetic fields around the Feshbach 
resonance. The behavior at 9.15 G is quite different from 
that at 9.05 G even though the numbers of atoms are 
almost the same. These results indicate that the scat- 
tering length ai2 significantly changes at B ev c± 9.1 G. 
The domain structures depicted in Fig.[2]are reproducible 
within the experimental fluctuations in the initial num- 
ber of atoms and the relative populations. 

We numerically solve coupled GP equations (1) for 
various values of Aa and if 12, and obtain the column 
densities of atomic clouds after the TOF (lower pan- 
els of Fig. [2j). The GP equations were solved by the 
Crank-Nicolson method. The numerically obtained col- 
umn densities are smoothed to account for the spatial 
resolution (3.2 /xm) of the CCD camera. The values of 



A a and Kyi are determined as follows. First, we assume 
Aa = and calculate the time evolution of the num- 
ber of atoms for various K12 using coupled GP equa- 
tions (1). We then compare the results with the exper- 
imental ones and estimate the value of Ky^. We con- 
firmed that the difference between the number of atoms 
for Aa = and that for Aa 7^ differ by less than 15 
%, which is comparable to the fluctuation in the initial 
number of atoms and does not affect the estimation of 
K12. We next calculate the column densities for var- 
ious values of Aa, compare them with the experimen- 
tal results, and estimate the value of Aa. For exam- 
ple, Fig. [3] shows a catalog of the column densities for 
K12 = 3 x 10 _13 cm 3 /s, which corresponds to B ev = 9.05 
G and 9.15 G in Fig. O Comparing the experimental 
column densities in Fig. [2] with numerical ones in Fig. [3j 
we can estimate scattering lengths. We assume that 
the most probable scattering l ength minimizes the root- 
mean-square deviation s = \J*YTz [ a exp(^) — <^cai(^)] 2 / n > 
where a exp (z) and a ca \(z) are optical densities in an ex- 
periment and a calculation at z, respectively. Figure 2] 
shows integrated optical densities obtained by the exper- 
iment at 9.05 G and numerical calculation at Aa = 7.5, 
8.0, and 8.5 ae- We find that the phase separation dy- 
namics is sensitive to the change in the scattering length 
by 0.50 ttB- The difference between experiments and cal- 
culations are expressed by s values as shown in Fig. [5] 
The value of Aa that minimizes the value of root-mean- 
square deviation s m i n is Aa = 8.0 an for B ev = 9.05 G. 
In the range of Aa > 8.0 ae, fluctuation of s is expected 
to be larger than that in Aa < 8.0 ae- This is because 
the density distribution of each component forms com- 
plex structures at Aa > 8.0 an while it shows simple 
domain structure in Aa < 8.0 an- The phase-separation 
dynamics can thus be used as a probe for estimating the 
scattering length in a binary BEC. 

The density patterns are sensitive to the differences 
in the scattering length for positive Aa. The accuracy 
of the estimated scattering lengths is therefore typically 
±1 ab for Aa > 3flB- However, in the close vicinity of 
the Feshbach resonance, the domain structure becomes 
moderate since the number of atoms decreases consid- 
erably, which makes accurate estimation of Aa difficult. 
For -7 ae ^ Aa < in Fig. [3j inhomogeneous density 
distribution is formed, even though the phase separation 
condition is not satisfied; we can estimate Aa = —6.5 an 
at 8 m i n for B ev = 9.15 G. This is because the number 
of 1 2) atoms decreases rapidly, giving rise to a nonequi- 
librium density distribution. Clear phase separation is 
necessary for precise estimation of the scattering length 
between |1) and |2). In our experiment, t ev < 25 ms is 
insufficient for observing the phase separation. For t ev > 
100 ms, the number of atoms in the |2) state becomes too 
small to estimate Aa. The evolution time of t ev = 75 ms 
is the most suitable to estimate the scattering length. 

For a uniform system, the Bogoliubov excitation spec- 
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FIG. 4: (Color online) Comparison between the experi- 
mentally observed density distributions at B ev = 9.05 G 
(dashed lines) and numerical ones (solid lines) obtained for 

(a) Aa = 7.5 aB, (b) 8.0 aB, and (c) 8.5 aB- 
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FIG. 5: (Color online) Root-mean-square deviation be- 
tween experiments and calculations at 9.05 G with K12 = 
3 x 10 -13 cm 3 /s. The dashed line shows 1.3 s m in, which de- 
fined as a threshold for determination of the error bar of the 
scattering length. 



trum has the form 3C 



(hu) 2 =e 



£+giini+g22n2±\(giin 1 -g22n2) 2 +4.n 1 n 2 gl2 



(4) 
where rij = \^j\ 2 is the atom density and e = ft k 2 /(2m) 
with k being the excitation wave number. When g 2 2 > 
911922, u is imaginary for e < [(gurit - ^ 2 2^>) 2 + 
4nin2gf 2 ] 1 ^ 2 — (gn n i +#22^2) and the system becomes 
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FIG. 6: (Color online) Estimated value of the change in the 
scattering length between the |1) and |2) states for a magnetic 
field around the Feshbach resonance. The scattering lengths 
at Smin are shown as filled squares, and fitted by Eq.Q (solid 
curve). The open squares correspond to s below the thresh- 
old value of 1.3 s m in- Theoretical prediction in Ref.[19] is 
indicated by the dashed curve. 



dynamically unstable against phase separation. The ex- 
ponential growth of the unstable mode is approximately 
given by 



exp 



/ dtlmuj{t) 



(5) 



where uj(t) depends on time since n\ and n^ decrease with 
time. The wavelength that maximizes Eq. (J5j) is most 
unstable, which is, for example, estimated to be 8.9 /xm 
for K\2 = 1 x 10 _13 cm 3 /s and Aa = 1.0 ae corresponding 
to 9.00 G, and 4.9 /im for K 12 = 3 x 10- 13 cm 3 /s and 
Aa = 8.0 <2b to 9.05 G in Fig. [2j respectively. These 
estimations of the most unstable wavelengths are in good 
agreement with the numerical solutions of Eq. (1) before 
the TOF. By the TOF expansion, the wavelengths in the 
density pattern become a few times larger. 

Figure [6] depicts the normalized scattering length Aa. 
The theoretical prediction in Ref. [19] is adjusted to the 
resonance magnetic field of 9.10 G obtained in our experi- 
ment (dashed curve). The scattering lengths are obtained 
by the method in FigO The filled squares correspond to 
Smin an d the open squares show accuracy ranges below 
1.3 8 m i n of the threshold. They are best fitted by Eq.Q 
(solid line) with B = 9.104 G (7b = 13 mG, AB = 3 
mG). The value of AB is in good agreement with that 
of the theoretical prediction and other experiments [19- 
l2ll [291 ] . The discrepancy in 7b may be caused by uncer- 
tainties near the resonance field. 

The accuracy of Aa in our method is comparable to 
that in the method using the Ramsey fringe [21]. At 
B ev > Bo with high density regime, the accuracy of 
Aa in our method is improved at large Aa because the 



phase separation dynamics is sensitive to Aa. On the 
other hand, the accuracy in the Ramsey fringe method 
at large Aa becomes worse owing to density inhomogene- 
ity caused by dramatic phase separation. We note that 
the discrepancies in the scattering length between ex- 
perimental estimation and theoretical prediction may be 
caused by entangled spin states or molecular states near 
the Feshbach resonance [lg, LL7[ • If these effects are taken 
into account, the accuracy of the method may be im- 
proved. 



V. CONCLUSIONS 

In conclusion, we observed the time evolution of bi- 
nary 87 Rb condensates in the |2, — 1) and |1, 1) hyperfine 
states around the Feshbach resonance at 9.10 G. In the 
vicinity of the Feshbach resonance, the miscibility of the 
two components is tuned to be both immiscible and mis- 
cible. Phase separation occurs on the lower field side of 
the Feshbach resonance, while miscible behavior is ob- 
served on the higher field side. We performed numerical 
simulations using coupled GP equations and proposed a 
new method for determination of the scattering length. 
We estimated the values of Aa and K12 by comparing 
the experimental and numerical density distributions of 
the atomic cloud. From systematic experiments and sim- 
ulations, we obtained the resonance curve of scattering 
length around the Feshbach resonance. 

Our method for determination of the scattering length 
can be used for not only spectroscopic states but also non- 
spectroscopic states, and is powerful technique for high 
density regime in both identical and different isotopes. 
Miscibility control via mixed-spin-channel Feshbach res- 
onance will open up new possibilities for multicomponent 
quantum fluids, such as controlled quantum phase transi- 
tion between miscible and immiscible phases with precise 
tuned scattering lengths. In addition, this technique can 
be applied for a precise measurement of a magnetic field 
below sub-milligauss range in the case when the scatter- 
ing length curve are well-known. 
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